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We present a measurement of the W boson mass in W — > ev decays using 1 fb _1 of data collected 
with the DO detector during Run II of the Fermilab Tevatron collider. With a sample of 499830 
W — > ev candidate events, we measure Mw = 80.401 ± 0.043 GeV. This is the most precise 
measurement from a single experiment. 

PACS numbers: 12.15.-y, 13.38.Be, 14.70.Fm 



Knowledge of the W boson mass (Mw) is currently a 
limiting factor in our ability to tighten the constraints on 
the mass of the Higgs boson as determined from internal 
consistency of the standard model (SM) [lj. Improving 
the measurement of Mw is an important contribution 
to our understanding of the electroweak (EW) interac- 
tion, and, potentially, of how the electroweak symmetry 
is broken. The current world-average measured value is 
M w = 80.399 ± 0.025 GeV Q from a combination of 
measurements from the ALEPH 0, DELPHI @|, L3 0], 



OPAL 3, DO @, and CDF 0,Q collaborations. 

In this Letter we present a measurement of Mw us- 
ing data collected from 2002 to 2006 with the DO detec- 
tor Q, corresponding to a total integrated luminosity of 
1 fb _1 (loj . We use the W — > ev decay mode because 
the DO calorimeter is well-suited for a precise measure- 
ment of electron energies, providing an energy resolu- 
tion of 3.6% for electrons with an energy of 50 GeV. 
The components of the initial state total momentum 
and of the neutrino momentum along the beam direc- 
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tion are immeasurable, so My/ is measured using three 
kinematic variables measured in the plane perpendicular 
to the beam direction: the transverse mass m T , the elec- 
tron transverse momentum pj,, and the neutrino trans- 
verse momentum p^- The transverse mass is defined as 
rriT = y / 2p? r pl^(l — cos A</>), where A</> is the opening 
angle between the electron and neutrino momenta in the 
plane transverse to the beam. The magnitude and direc- 
tion of Pt are inferred from the event missing transverse 
energy $ T ). The My/ measurement is made by compar- 
ing data spectra of tut, p^i and^ T with probability den- 
sity functions (templates) for these spectra constructed 
from Monte Carlo simulation with varying input My/ 
values. 

The DO detector [§| contains tracking, calorimeter, 
and muon systems. Silicon microstrip tracking detectors 
(SMT) near the interaction point cover pseudorapidity 
\r)\ < 3 to provide tracking and vertex information. The 
central fiber tracker surrounds the SMT, providing cov- 
erage to | X] | « 2. A 2 T solenoid surrounds these track- 
ing detectors. Three uranium, liquid-argon calorimeters 
measure particle energies. The central calorimeter (CC) 
covers \rj\ < 1.1, and two end calorimeters (EC) extend 
coverage to \r]\ « 4. The CC is segmented in depth 
into eight layers. The first four layers are used primarily 
to measure the energy of photons and electrons and are 
collectively called the electromagnetic (EM) calorimeter. 
The remaining four layers, along with the first four, are 
used to measure the energy of hadrons. Intercryostat 
detectors (ICD) provide added sampling in the region 
1.1 < \r]\ < 1.4 where the CC and EC cryostat walls 
degrade the calorimeter energy resolution. A three level 
trigger system selects events for recording with a rate of 
100 Hz. 

Events are initially selected using a trigger requiring 
at least one EM cluster found in the CC with transverse 
energy threshold varying from 20 GeV to 25 GeV de- 
pending on run conditions. Additionally, the position of 
the reconstructed production point of a W or Z boson 
along the beam line is required to be within 60 cm of the 
center of the detector. 

Candidate W boson events are required to have one 
EM cluster reconstructed in the CC, with pj, > 25 GeV 
and |?7| < 1.05 where r\ is the pseudorapidity measured 
with respect to the center of the detector. The EM clus- 
ter must pass electron shower shape and energy isola- 
tion requirements in the calorimeter, be within the cen- 
tral 80% of the electromagnetic section of each CC mod- 
ule, and have one track matching in (77, </>) space, where 
the track has at least one SMT hit and pr > 10 GeV. 
The central 80% requirement is applied to the <fi coor- 
dinate only and excludes regions with slightly degraded 
energy resolution. The event must satisfy Ifi T > 25 GeV, 
u T < 15 GeV, and 50 < m T < 200 GeV. Here $ T is 
the magnitude of the vector sum of the transverse en- 



ergy of calorimeter cells above read out threshold, exclud- 
ing those in the coarse hadronic layer and in the inter- 
cryostat detector, and ut is the magnitude of the vector 
sum of the transverse component of the energies mea- 
sured in calorimeter cells excluding those associated with 
the reconstructed electron. This selection yields 499,830 
candidate W — > ev events. Throughout this Letter we 
use "electron" to imply either electron or positron. 

We use Z — > ee events for calibration. Candidate Z bo- 
son events are required to have two EM clusters satisfying 
the requirements above. Both electrons must have p^, > 
25 GeV. One must be reconstructed in the CC and the 
other in either the CC or EC (1.5 < < 2.5). The as- 
sociated tracks must be of opposite charge. Events must 
also have u T < 15 GeV and 70 GeV < m ee < 110 GeV, 
where m ee is the invariant mass of the dielectron pair. 
Events with both electrons in the CC are used to de- 
termine the EM calibration. There are 18,725 candidate 
Z — » ee events in this category. 

The backgrounds in the W boson sample are Z — * ee 
events in which one electron escapes detection, multi- 
jet events (MJ) in which a jet is misidentified as an 
electron with Jp T arising from misreconstruction, and 
W — ► tv — > evvv events. The background from Z bo- 
son events arises from electrons which traverse the gap 
between the CC and EC. The tracking efficiency in this 
region is high, so this background is estimated by select- 
ing data events passing the W boson selection in which an 
additional track is pointing at the gap region. The MJ 
background is determined using a sample obtained by 
removing the track matching requirement for the elec- 
tron candidates. The probabilities for background and 
W boson signal events in this sample to have a match- 
ing track are measured in control samples. The number 
of events in the sample without the track requirement 
and the two probabilities are then used to determine the 
number of MJ background events in the final W boson 
sample. The W — > tv — > evvv contribution is deter- 
mined from detailed simulation of the process using the 
DO GEANT [ijj-based simulation. The backgrounds ex- 
pressed as a fraction of the final sample are (0.90±0.01)% 
from Z -> ee, (1.49±0.03)% from MJ, and (1.60±0.02)% 
from W — > tv — > evvv. 

W and Z boson production and decay kinematics are 
simulated using the RESBOS [ijj next-to-leading order 
generator which includes non-perturbative effects at low 
boson pt- These effects are parametrized by three con- 
stants (gi, q-y and 53) whose values are taken from global 
fits to data [13| . The radiation of one or two photons is 
performed using the PHOTOS [3| program. 

Detector efficiencies and energy response and resolu- 
tion for the electron and hadronic energy are applied 
to the resbos+photos events using a fast paramet- 
ric Monte Carlo simulation (fastmc) developed for this 
analysis. The fastmc parameters are determined using 
a combination of detailed simulation and control data 



5 



samples. The primary control sample used for both the 
electromagnetic and hadronic response tuning is Z — > ee 
events. W boson events are also used in a limited man- 
ner, as are events recorded in random beam crossings, 
with or without requiring hits in the luminosity coun- 
ters. 

Since the Z boson mass and width are known with high 
precision from measurements [IH at the CERN e + e~ 
collider (LEP), these values are used to calibrate the 
electromagnetic calorimeter response assuming a form 
^moas _ a ^truo _j_ p with a and P constants determined 
by calibration. The Mw measurement presented here is 
effectively a measurement of the ratio of W and Z boson 
masses. Figure Q] shows a comparison of the m ee distribu- 
tions for data and fastmc, as well as the x distribution 
defined as the difference between data and the fastmc 
prediction divided by the statistical uncertainty on the 
difference. 

The other major calibration is that of the hadronic en- 
ergy in the event, which includes energy recoiling against 
the boson. The hadronic response (resolution) is tuned 
using the mean (width) of the ryi m b distribution in Z — > ee 
events in bins of p|f . Here ryimb is defined as the sum of 
the projections of the dielectron momentum (p!p ) and ut 
vectors in the transverse plane on the axis bisecting the 
dielectron opening angle [la ]. 




:T i I 1 1 

70 75 80 85 90 95 100 105 110 
m ee (GeV) 



FIG. 1: (a) The dielectron invariant mass distribution in Z — *■ 
ee data and from the fast simulation fastmc and (b) the 
X values where Xi — [N% — (FASTMCi )]/(?% for each point in 
the distribution, Ni is the data yield in bin i and ai is the 
statistical uncertainty in bin i. 

To determine Mw, fastmc template distributions for 
rriT, Vt-> and-^y are generated at a series of test Mw val- 
ues at intervals of 10 MeV with the backgrounds added 
to the simulated distributions. A binned likelihood be- 
tween the data and each template is then computed. The 
resulting log likelihoods as a function of mass arc fit to 
a parabola. The minimum point of the parabola defines 
the measured Mw value. The fits are performed sepa- 
rately for each of the tot, Vti and^ T distributions, and 
the fit ranges were chosen to minimize the total expected 
uncertainty on for each distribution. 



A test of the analysis procedure is performed using 
events produced by the detailed GEANT Monte Carlo sim- 
ulation treated as collider data. The methods used for 
the data analysis are applied to the simulated events, in- 
cluding the fastmc tuning using the simulated Z — > ee 
events. Each of the Mw fit results using the mr, Pj>> an d 
Ifirp distributions agree with the input Mw value within 
the 20 MeV total uncertainty of the test arising from 
Monte Carlo statistics. 

During the fastmc tuning performed to describe the 
collider data, the Mw values returned from fits are 
blinded by the addition of an unknown constant offset. 
The same offset was used for %, pfp and tfJ T . This 
allowed the full tuning on the W and Z boson events 
and internal consistency checks to be performed without 
knowledge of the final result. Once the important data 
and fastmc comparison plots have acceptable x distribu- 
tions, the results are unblinded. The Z boson mass value 
from the post-tuning fit is 91.185 ± 0.033 (stat) GeV, in 
agreement with the world average of 91.188 GeV used for 
the tuning. The Mw results from data after unblinding 
are given in Table [IJ The tyit , p^ , and lfi T distributions 
showing the data and fastmc template with background 
for the best fit Mw are shown in Fig. O 



TABLE I: Results from the fits to data. The uncertainty is 
only the statistical component. 



Variable Fit Range (GeV) 


M w (GeV) 


x7dof 


mr 65 < mr < 90 
Pt 32 < p e T < 48 
$ T 32<# T <48 


80.401 ± 0.023 
80.400 ± 0.027 

80.402 ± 0.023 


48/49 
39/31 
32/31 



The systematic uncertainties in the Mw measurement 
arise from a variety of sources, and can be categorized as 
those from experimental sources and those from uncer- 
tainties in the production mechanism. The systematic 
uncertainties are summarized in Table HT1 

The uncertainties on the electron energy calibration 
and the hadronic recoil model are determined by simul- 
taneously varying the parameters determined in the tun- 
ing to Z — > ee events by one statistical standard de- 
viation including correlation coefficients. The electron 
energy resolution systematic uncertainty is determined 
by varying resolution parameters determined in the fit to 
the width of the observed Z — > ee m ee distribution. The 
shower modeling systematic uncertainties are determined 
by varying the amount of material representing the de- 
tector in the detailed simulation within the uncertainties 
found by comparing the electron showers in the simula- 
tion to those observed in data. No effect was seen when 
studying possible systematic bias for the energy loss dif- 
ferences arising from the differing E or rj distributions for 
the electrons from W and Z boson decay. The quoted 
systematic uncertainty is due to the finite statistics of the 
event samples from the tuned detailed simulation that are 
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FIG. 2: The (a) tut, (b) py, and (c) Ifi T distributions for data and FASTMC simulation with backgrounds. The \ values are 
shown below each distribution where \i — [Ni — (fastmc^)]/ cti for each point in the distribution, Ni is the data yield in bin i 
and only the statistical uncertainty is used. The fit ranges are indicated by the double-ended horizontal arrows. 



TABLE II: Systematic uncertainties of the Mw measurement. 

AM W (MeV) 



Source 




Pt 




Electron energy calibration 


34 


34 


34 


Electron resolution model 


2 


2 


3 


Electron shower modeling 


4 


6 


7 


Electron energy loss model 


4 


4 


4 


Hadronic recoil model 


6 


12 


20 


Electron efficiencies 


5 


6 


5 


Backgrounds 


2 


5 


4 


Experimental Subtotal 


35 


37 


41 


PDF 


10 


11 


11 


QED 


7 


7 


9 


Boson pr 


2 


5 


2 


Production Subtotal 


12 


14 


14 


Total 


37 


40 


43 



divided into statistically independent categories based on 
instantaneous luminosity, time, the total hadronic trans- 
verse energy in the event, the vector sum of the hadronic 
energy, and electron pseudorapidity range. The fit ranges 
are also varied. The results are stable to within the mea- 
surement uncertainty for each of these tests. 

The results from the three methods have combined 
statistical and systematic correlation coefficients of 0.83, 

0.82, and 0.68 for (rnr, Pt)i ( to t,-^t)' ano - (Pt>-^t) re ~ 
spectively. The correlation coefficients are determined 
using ensembles of simulated events. The results are com- 
bined [2l| including these correlations to give the final 
result 

M w = 80.401 ± 0.021 (stat) ± 0.038 (syst) GeV 
= 80.401 ± 0.043 GeV. 



used to transport calibrations from the Z to the W sam- 
ple. The electron efficiency systematic is determined by 
varying the efficiency by one standard deviation. Table ITT1 
also shows the Mw uncertainties arising from variation 
of the background uncertainties indicated above. 

Among the production uncertainties, the parton dis- 
tribution function (PDF) uncertainty is determined by 
generating W boson events with the PYTHIA 17| pro- 
gram using the CTEQ6.1M [lgj] PDF set. The CTEQ 
prescription [l8| is used to determine a one standard de- 
viation uncertainty Q on Mw. The QED uncertainty is 
determined using WGRAD [lj| and ZGRAD [2(J, varying 
the photon-related parameters and assessing the varia- 
tion in Mw and by comparisons between these and PHO- 
TOS. The boson px uncertainty is determined by varying 
gi by its quoted uncertainty [13] . Variation of g\ and g% 
has negligible impact. 

The quality of the simulation is indicated by the good 
X 2 values computed for the difference between the data 
and FASTMC shown in the figures. The data are also sub- 



The dominant uncertainties arise from the available 
statistics of the W — > ev and Z — ► ee samples. Thus, 
this measurement can still be expected to improve as 
more data are analyzed. The Mf measurement reported 
here agrees with the world average and the individual 
measurements and is more precise than any other single 
measurement. Its introduction in global electroweak fits 
is expected to lower the upper bound on the SM Higgs 
mass, although it is not expected to change the best fit 
value [l[. 
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